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ANALYSIS OF THE GAMMA-RAY- BACKSCATTER 
TECHNIQUE FOR DIRECT MEASUREMENT OF THE DENSITY 
OF THE MARTIAN ATMOSPHERE 

By Peter J. LeBel 
Langley Research Center 

SUMMARY 

Studies have been conducted of a technique that uses the backscattering of gamma 
radiation emitted by a radioisotope source to measure the ambient density of the Martian 
atmosphere during entry of a spacecraft into the atmosphere. The results show the fea- 
sibility of the measurement up to an altitude of about 100 000 feet (30.5 km) in the least 
dense of the Voyager-Mars (VM) model atmospheres. No probes or special windows 
through the spacecraft skin or heat shield are required because of the penetrating ability 
of the radiation. A source strength at Mars of 20 curies (l curie = 3.7 x lO^O disintegra- 
tions per second) of gamma radiation is required. The measurement is essentially unaf- 
fected by the composition of the atmosphere, by the shock layer, and by motions of the 
vehicle. An engineering- model density sensor, which has been fabricated and tested, has 
a mass of 5 pounds (2.3 kg), a volume of 90 inches3 (1475 cm3), and requires 5 watts of 
power for operation. Further miniaturization of the device appears practical. Steriliza- 
tion of the sensor is not expected to present major problems. 

INTRODUCTION 

Density is a quantity of prime importance to be measured during the probing of plan- 
etary atmospheres. The atmospheric density profile is not only of considerable scientific 
value but is also an important engineering measurement needed for optimum aerodynamic 
design of future spacecraft. In general, information on the atmospheric density from a 
planet's surface to as high in its atmosphere as possible is desired. However, the mission 
profile usually presents limitations. For a large part of the proposed Martian atmosphere 
entry trajectories, the instrument package will be located behind an ablative heat shield 
and the entry capsule engulfed in a high-temperature flow field. Under these conditions, 
measurement techniques which avoid ports, probes, or special thin windows that might 
compromise the integrity of the heat shield or structure are desirable. In addition, the 
shock layer and vehicle motions should have a minimum effect on the measurement. 

A technique employing backscattering of gamma radiation emitted by a radioisotope 
source to measure density has been shown to be effective under these conditions. This 



concept is shown schematically in figure 1. It consists of a source of gamma radiation, 
a detector, and shielding material to eliminate direct transmission between the two. 
Because of their penetration ability, gamma rays can propagate out through the vehicle 
structure and heat shield without severe attenuation. A fraction of the emerging gamma- 
ray flux is backscattered by the atmosphere into the detector. The backscatter interac- 
tion, Compton scattering of gamma rays by atomic electrons, is linearly dependent upon 
the atmospheric density in the scattering region. By proper location of sensor compo- 
nents and collimation of the source emission and detector field of view, the effective 
scattering volume or measurement region can be located in the undisturbed ambient 
gas beyond the shock wave. Calibration accounts for any attenuation of the gamma-ray 
beam in the heat shield; attenuation in the shock layer is generally very small and can 
be neglected. The density profile is then obtained by correlating the measured density 
with spacecraft altitude during entry. 

For measurement of Martian atmospheric density, the potential advantages of the 
gamma-ray-backscatter technique are 

(1) The measurement is a direct, linear function of the ambient atmospheric 
density 

(2) No probes or special windows through the vehicle skin or heat shield are 
required because of the penetrating ability of the radiation 

(3) Atmospheric composition measurements are not required for accurate 
density measurements 

(4) The shock layer and vehicle motions have a negligible effect upon the 
measurement. 

In general, the gamma-ray-backscatter technique is applicable to other planetary 
atmospheres; however, this report will discuss the results of NASA in-house and con- 
tractor studies of the feasibility of using the technique to measure the ambient density 
of the Martian atmosphere during entry. 

SYMBOLS 

A effective atomic weight of the atmosphere 

B background counting rate 

Bp) direct-transmission background counting rate 
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Bjg natural-radiation background counting rate 

Bg skin-scatter background counting rate 

f 1 ,f 2 3 ,f4 factors of K which express differences in the various sensor-response 
parameters between the time of sensor calibration and the time of mea- 
surement during atmospheric entry (See appendix.) 

G gross or total counting rate 

K proportionality constant relating backscattered counting rate and atmospheric 

density 

N Avogadro's number 

R counting rate per unit density measured during sensor calibration 

S backscattered counting rate, or signal 

tg time for background measurement 

tQ time for gross or total counting-rate measurement 



°s 


a 


s,e 


effective atomic number of the atmosphere 
density of the scattering medium (atmosphere) 
standard deviation 

total atmospheric Compton scattering cross section 

Compton scattering cross section per electron (a function of the gamma-ray 
energy) 
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BACKGROUND 


Gamma-Ray Scattering 

A discussion of the several ways in which a gamma ray can interact with matter 
can be found in the literature. (See, for example, ref. 1.) Only four interactions are 
significant for gamma-ray energies in the range between 10 keV and 10 MeV. These 
interactions are illustrated in figure 2, which is a plot of the relative probability of 
occurrence of each interaction as a function of gamma-ray energy, for nitrogen (N2) or 
carbon dioxide (CO2) at standard temperature and pressure. (Curves for N2 are pre- 
sented in ref. 2 from cross-section data tabulated in ref. 3. Similar calculations for 
CO2 and other gases from ref. 3 data have been found to yield similar results.) The 
photoelectric effect is an ionization process in which the gamma- ray energy is completely 
transferred to an electron bound to the nucleus of an atom. In the same general energy 
region is Rayleigh, or coherent, scattering. This interaction, involves small-angle scat- 
tering of the gamma ray from tightly bound atomic electrons and has a low probability of 
occurrence. At higher energies pair production becomes important. In this phenomenon, 
the gamma ray interacts with the intense field surrounding the nucleus and is converted 
into a positron-electron pair. A minimum gamma-ray energy of 1.02 MeV, twice the 
rest-mass energy of an electron, is required. Figure 2 also illustrates Compton, or 
incoherent scattering, in which the gamma ray acts like a particle and undergoes a col- 
lison with an electron. The electron absorbs some of the gamma-ray energy and is 
freed from the atom. The scattered gamma-ray has an energy which depends upon its 
original energy, the energy absorbed by the electron and the scattering angle. As can 
be seen in figure 2, over much of the energy range shown, Compton interactions have 
the highest probability of occurrence. Also, of the four mechanisms involved, only 
Compton scattering can produce a large-angle or backscattered gamma ray which will 
reach the detector. 


Concept- Feasibility Flight Tests 

Two flight tests have been conducted in the earth's atmosphere to examine the fea- 
sibility of using backscattered gamma rays to measure atmospheric density. The launch 
vehicle for these tests was the Nike-Apache sounding rocket. The sensor payload, shown 
schematically in figure 3, included a source of approximately 30 curies of the isotope 
cerium-144 (l curie = 3.7 x 10*0 disintegrations per second) and a sodium iodide scintil- 
lation detector. 

Test data for one of the flights are shown in figure 4. The counting rate and den- 
sity are presented as a function of altitude. For comparison, density determined from 
Arcasonde measurements at low altitudes as well as density data from reference 4 at 
higher altitudes are shown. The point scatter at high altitudes occurs because of the 
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statistical uncertainty inherent in counting random nuclear events, and also because a 
large background counting rate had to be subtracted from each data point. 

The flight-test results demonstrated the feasibility of using backscattered gamma 
rays to measure atmospheric density but indicated that for measurement of a thin atmo- 
sphere such as that of Mars, improvements in signal-to-noise ratio were necessary. As 
a result, investigations were initiated to reduce background noise and to develop a high- 
purity radioisotope source. A study of the sterilization resistance of a miniaturized 
gamma- ray -backscatter density sensor was also initiated. Additional details of the fea- 
sibility flight tests can be found in reference 5. 

RESPONSE OF MARTIAN ATMOSPHERE SENSOR 


Signal Response 

In studying the feasibility of using backscattered gamma rays to measure the density 
of the Martian atmosphere, the parameters which affect sensor response must be exam- 
ined. For a typical measurement geometry such as that shown in figure 1, the backscat- 
tered counting rate, or signal, S is a function of the gamma-ray flux emerging from the 
source, gamma-ray energy, attenuation in the gamma- ray path, Compton scattering cross 
section in the atmosphere, detector area and efficiency, and the geometry parameters, 
such as the source-detector separation and the collimation angles. For fixed source, 
detector, and geometry, all these parameters are constant with two exceptions, the atten- 
uation in the shock layer and the atmospheric Compton scattering cross section 6 S . As 
discussed in the appendix, for the Martian atmosphere the attenuation in the shock layer, 
although not constant, is expected to be very small and can be neglected. The total atmo- 
spheric Compton scattering cross section is given by 


a s ~ °s,e 



where 6 s>e and N are constant. 

be represented by 

* 


Thus, the backscattered counting rate, or signal, can 
S = K P | 


where K is a proportionality constant which includes the source, detector, and geometry 
parameters previously mentioned. 

This expression shows that the response of the device is a direct function of atmo- 
spheric density. The only other atmospheric parameter which affects the signal is the 
ratio of atomic number to atomic weight Z/A. This ratio is equal to 0.5 for all gases in 
the Voyager-Mars model atmospheres, VM-1 through VM-10 (table 1 from ref. 6) with the 
exception of argon which has a Z/A ratio of approximately 0.45. This value represents 
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a 3.5 percent uncertainty in density measurement for the maximum argon concentration 
postulated (VM-6 model atmosphere). Thus, atmospheric composition measurements are 
not required for accurate density measurements. 

Background Noise 

During descent of the sensor through the Martian atmosphere, background radiation 
will contribute to the total counting rate. It is important to reduce background noise for 
maximum measurement sensitivity and to precisely measure background counting rate 
for maximum measurement accuracy. The background counting rate will be measured 
just prior to atmospheric entry and subtracted from the total counting rate so that the 
backscattered counting rate, or signal, can be determined. Since the signal increases 
with increasing density and the background noise will be approximately constant, the 
signal-to- noise ratio will continually improve during entry. Three factors which con- 
tribute to the background noise are 

(1) Natural- background radiation. In reference 7 it was estimated that cosmic rays 
and other sources of natural-background radiation would contribute on the order of 

50 counts per second to the total counting rate for a typical detection system and mea- 
surement geometry. 

(2) Direct-transmission noise. Direct transmission of gamma rays between source 
and detector can be reduced by shielding; however, high-energy radiation is difficult to 
shield. In the flight payloads described earlier, a thickness of 15 inches (38 cm) of lead 
was required to shield the detector from a small amount of 2-MeV radiation from the 
cerium- 144 source. If the source had only emitted the primary energy, 134 keV, less than 
0.4 inch (1 cm) of shielding would have been required. Since weight is of prime impor- 
tance on a Mars mission, the radiation source must have no high-energy gamma-ray com- 
ponent which would require extensive shielding. Additional discussion of radiation source 
requirements can be found later in this paper. 

(3) Skin-scatter noise. The background counting rate which was measured on the 
Nike-Apache flight tests was about an order of magnitude greater than had been predicted 
for cosmic radiation and other natural sources. As discussed in reference 8, the major 
part of this high background noise level was found to be caused by multiple scattering of 
gamma rays within the vehicle skin and into the detector. Results of in-house studies and 
those reported in reference 8 show several potential methods of minimizing skin-scatter 
noise. It can be reduced considerably by choosing a geometry that requires the gamma 
ray to undergo a third scatter before being detected. Skin-scatter noise has also been 
found to be a strong function of the source-detector separation distance. Small increases 
in this separation will significantly reduce the noise without appreciably affecting the back- 
scattered signal. (See ref. 9.) 
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Measurement Range 


For a typical geometry, such as that shown in figure 1, sensor response to atmo- 
spheric density can be estimated by examining the probability functions that determine 
the path of any gamma ray emitted by the source, scattered by the atmosphere toward the 
detector, and detected. The product of these probabilities integrated over the total volume 
of interaction represents the response of the device. This approach is further discussed 
in reference 8. 

For a typical measurement of Martian atmospheric density, sensor response was 
estimated by using the techniques of reference 8, the geometry of figure 1, and the fol- 
lowing conditions: 

(1) Source-detector separation, 1.5 feet (0.5 m) 

(2) Source energy, 100 keV 

(3) Attenuation of gamma- ray beam in skin, 10 percent; this is what might be 
expected in a low-density silicone-elastomer heat-shield material of the type being con- 
sidered for Mars entry missions 

(4) Detector, 3-inch (7.6-cm) diameter, 100-percent efficient for 100-keV gamma 
radiation 

(5) Density varied with altitude for the extremes of Mars model atmospheres, VM-8 
and VM-9, as shown in figure 5 

(6) Source strength at Mars of 20 gamma curies (a gamma curie is equivalent to 
3.7 X 10l0 gamma rays per second) 

(7) The capsule heat shield was discarded at an altitude of 25 000 feet (7.6 km) for 
both VM-8 and VM-9 

For these conditions, sensor response in counts per second was found to be 

S = (7.0 X 1 0 8 )p 


before heat- shield ejection, and 

S = (7.6 x 10 8 )p 

after heat- shield ejection, where p is in units of g/cm 8 . 

By using the density- altitude data of figure 5 and 1- second counting intervals, sen- 
sor response is shown as a function of altitude for VM-8 and VM-9 in figure 6. For VM-8, 
the device begins to respond to density at an altitude of about 150 000 feet (45.8 km) with 
a counting rate of a few counts per second increasing to 10 000 counts per second at the 
Martian surface. In the much more dense VM-9 atmosphere, the sensor begins to respond 
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to density at an altitude of 400 000 feet (122 km) with counting rate increasing to 
20 000 counts per second at the planet's surface. 

If necessary, to avoid interference with the landed spacecraft experiments, the 
gamma- ray backscatter sensor could be ejected from the spacecraft when the heat shield 
is discarded, or a shielding mechanism could be actuated to "turn off" the source. 


Measurement Accuracy 

In the appendix, an expression is derived for the relative standard deviation or 
uncertainty in density in terms of the total or gross counting rate G, background counting 
rate B, net counting rate, or signal, S, and proportionality constant K as 


g(p) 

P 


cr 2 (G) ct 2 (B) a 2 (K) /g 2 + B 2 \ 
S 2 S 2 K 2 \ S 2 / 


The various factors which influence the standard deviations of G, B, and K are also 
discussed in the appendix. 

The precision of a measurement of Martian atmospheric density can be estimated 
by using this expression and the following assumptions: 

(1) Signal varies with altitude as shown in figure 6 and discussed in the preceding 
section 

(2) The total background counting rate is 100 counts per second. Fifty counts per 
second are due to skin scattering and the remainder comes from such sources as natural 
background and direct transmission. The background reduces to 50 counts per second 
after heat- shield ejection, as shown in figure 6. 

(3) The systematic error is 0.03 (3 percent). A quantitative analysis of all 

K. 

factors which affect is beyond the scope of this report. However, based on a study 

K 

of the feasibility flight-test results, an estimate of 3 percent systematic error appears 
reasonable. 

The resultant 1-sigma uncertainty in density measurement is listed in table 2 and 
plotted in figure 7. The results show for VM-8 a 100 percent uncertainty at 135 000 feet 
(41.2 km), 10 percent at 90 000 feet (27.5 km), and 3 percent at the Martian surface. For 
VM-9 comparable results are 100 percent at 370 000 feet (112.8 km), 10 percent at 
235 000 feet (71.7 km), and 3 percent at the surface. 
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SYSTEM REQUIREMENTS 


Source 

A primary characteristic of the gamma- ray source is the energy spectrum of the 
emitted gamma rays. Energy selection for a backscatter sensor is complicated by having 
several conflicting criteria. As described in reference 2, these criteria include aware- 
ness that in order to use the technique, the gamma rays must penetrate the heat shield and 
return after having lost some energy in scattering. Maximum penetration ability implies 
the need for a high energy. However, the cross section for scattering at large angles is 
a maximum at low energy. A low gamma- ray energy is also desirable so that weight 
required to shield the detector from the source does not become impractical as previously 
discussed. Evaluation of these conflicting criteria can best be accomplished by a trade- 
off exercise. For a measurement of the Martian atmosphere, an optimum range of 
source energy is 50 keV to about 150 keV. (See ref. 10.) 

The most ideal source would be monoenergetic; thus all the gamma rays would be 
emitted at the energy of interest, and a narrow energy discrimination window could be 
established to increase the signal-to-noise ratio. The radioisotope comes closest to 
being the ideal monoenergetic source. The inherent high reliability and low weight com- 
pared with an X-ray tube make the radioisotope most suitable for the Mars application. 

Desirable source characteristics for the measurement of Martian atmospheric den- 
sity are summarized as follows: 

(1) Gamma- ray- emitting radioisotope 

(2) Primary gamma- ray energy between 50 and 150 keV 

(3) No high-energy gamma ray or beta particle which would require excessive 
shielding mass to minimize direct transmission noise 

(4) Source half-life, minimum of 150 to 200 days 

(5) High specific activity (curies per gram of source material) 

Under interagency agreement with Langley, the Oak Ridge National Laboratory of 
the Atomic Energy Commission has identified gadolinium- 153 as the most suitable radio- 

I cn 

isotope source for the Martian atmospheric measurement. (See ref. 11.) 10J Gd has a 
half-life of 241 days and emits gamma rays with an energy of about 100 keV and an abun- 
dance of 46 percent. It can be produced by two methods. The first is to irradiate a tar- 
get enriched with j n a nuc i ear reactor. l^Gd i S produced by the reaction 

152 Gd + neutron ► 153 Gd + gamma ray 


■M 
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Material produced by this method has a low specific activity and high cost because of the 
requirement for a highly enriched target. Also, terbium impurities (-^Tb emits high- 
energy gamma rays) were found in irradiated samples of gadolinium. Terbium- gadolinium 
chemical separations are difficult, and it is doubtful that a high-purity source could be 
fabricated. 

The second technique is to irradiate a natural europium target, containing about 
50 percent ^^Eu, in a reactor. The reaction sequence for production of *^Gd is 

■^Eu + neutron ► gamma ray + *®^Eu j 3e ^ a dec ay ^ ^Gd + neutron *- *^Gd + gamma ray 

The isotope produced by this method has a high specific activity and can be readily puri- 
fied to less than 10 parts per million of contamination. Also, because the target for irra- 
diation does not require enrichment, the isotope can be produced fairly inexpensively. 

The production sequence for obtaining the desired quantities of gadolinium is now being 
optimized. (See ref. 11.) 

For a specific mission, the required source strength depends upon a number of 
parameters, including the desired source output at time of use, source self-absorption, 
source half-life, and abundance of the gamma ray of interest. For example, to obtain an 
output of 20 gamma curies during descent through the Martian atmosphere, it would be 
necessary to initially produce 150 to 200 curies of ^^Gd. 

Detection System 

In order to minimize the required source strength, the detection system must be 
designed for maximum gamma detection efficiency. This requirement dictates the need 
for a scintillation-crystal — photomultiplier-tube detector. The detector area should 
be maximized consistent with such factors as the increased shielding weight and the 
increased susceptibility of large crystals to physical damage from shock and vibration. 

The detector output can be processed by using fairly standard techniques to produce 
a signal compatible with the telemetry system. Either analog or digital techniques can 
be used. 

Regardless of the processing network employed, some form of automatic gain con- 
trol (AGC) will be necessary to compensate for the poor gain stability of a scintillation 
detector with temperature. A standard AGC technique in nuclear systems is to use the 
Counting rate derived from a small, long half-life gamma or alpha radiation source placed 
near the crystal to control the gain by changing the photomultiplier tube voltage. 

Physical Characteristics 

Based on conceptual design studies (ref. 7) updated to take into account the 
latest knowledge of the Martian atmosphere and recent Mars mission studies, an 
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engineering-model density sensor has been fabricated for the purpose of sterilization 
testing. The device has a mass of 5 pounds (2.3 kg), a volume of 90 inches^ (1475 cm3), 
and requires 5 watts of power for operation. These figures include the radioisotope 
source, all necessary shielding, scintillation detector, signal conditioning electronics, 
and mounting hardware. Further miniaturization of the device appears practical. 

STERILIZATION 

The requirement for the sensor to withstand and be unaffected by sterilization pre- 
sents constraints in the selection of materials and components. Studies have shown that 
components of the type required can be found which have been qualified in the environment 
specified by reference 12. The major area of uncertainty has been the availability of 
suitable photomultiplier tubes. However, development programs being conducted by pho- 
totube manufacturers indicate that suitable devices can be obtained. (See, for example, 
ref. 13.) Testing of the engineering-model density sensor mentioned above is planned to 
verify the resistance of the device to the sterilization environment. 

CALIBRATION 

Calibration is necessary to establish sensor response as a function of atmospheric 
density. This must be accomplished in a vacuum chamber where gas density can be care- 
fully controlled and varied. Pressure and temperature measurements are made to moni- 
tor gas density, and counting- rate measurements are made to monitor sensor output. 

The calibration should take place in as large a vacuum chamber as possible to 
minimize scattering of gamma rays from the chamber walls and into the detector. In 
the feasibility flight-test program, calibration tests were conducted in the largest avail- 
able chamber, the 60-foot vacuum sphere at the Langley Research Center. Even so, it 
was necessary to perform the tests with an 0.5-curie source, since wall scattering would 
have exceeded the counting-rate capability of the sensor with the 30-curie flight source. 
The resultant calibration data were then scaled by the ratio of the source strengths, 

30 curies/0.5 curie. 

As discussed in the appendix, this scaling technique requires that the ratio of the 
gamma- ray outputs of the two sources be precisely known. To accomplish this, both 
calibration and flight sources are designed to have the same self-absorption and contain 
radioactive material from the same production batch. Then, to verify that the sources 
are constructed identically, direct-transmission measurements of the gamma- ray output 
of each source are made. 
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RADIATION SAFETY 


Handling and use of radioisotope sources of the type and size required for this sen- 
sor requires certain precautions to satisfy applicable regulations and minimize radiation 
exposure to personnel. Acceptable procedures for the safe handling of a large gamma- 
ray source developed and exercised during the feasibility flight-test program are directly 
applicable to a planetary mission. In addition, the engineering-model sensor mentioned 
earlier was designed with a remotely actuated shield mechanism to cover the source 
during prelaunch operations so that the radiation dose to individuals working in the vicin- 
ity of the sensor is below the limit permitted in an unrestricted area of 2 millirem per 
hour, 0.5 rem per year. 

The source-shield assembly will also incorporate sufficient thermal protection to 
protect the structural integrity of the source in case of a fireball explosion during launch 
operations. 


CONCLUDING REMARKS 

Studies have been conducted of the technique employing backscattering of gamma 
radiation emitted from a radioisotope source to directly measure the ambient density of 
the Martian atmosphere during entry of a spacecraft into the atmosphere. The measure- 
ment is essentially unaffected by the composition of the atmosphere, by the shock layer, 
and by motions of the vehicle. The results show the feasibility of the measurement during 
high-speed entry through the heat shield as well as after heat-shield ejection. For a typ- 
ical configuration, the estimated 1-sigma uncertainty in density measurement is ±15 per- 
cent at 100 000 feet (30.5 km) altitude in a VM-8 (least dense) model atmosphere, 
improving to ±3 percent at the Martian surface. These estimates were made for a 
source strength of 20 curies (l curie = 3.7 x 10*0 disintegrations per second) of gamma 
radiation at time of measurement. This would require production of 150 to 200 curies of 
gadolinium- 153, the most suitable source for the Martian atmospheric measurement. 
Techniques have been developed for producing this isotope with the proper purity. An 
engineering-model density sensor, which has been fabricated and tested, has a mass of 
5 pounds (2.3 kg), has a volume of 90 inches3 (1475 cm3), and requires 5 watts of power 
for operation. Further miniaturization of the device appears practical. Tests are 
planned to verify the sterilization resistance of the engineering-model sensor. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 7, 1969. 
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APPENDIX 


DERIVATION OF DENSITY- ERROR EXPRESSION 


An expression for the uncertainty in density measurement as a function of total and 
background counting rates and of systematic error is derived. 

By assuming that the backscattered counting rate, or signal, is directly proportional 
to atmospheric density, S = Kp or 

P = | (Al) 

The signal is obtained by subtracting the background counting rate from the total or 
gross counting rate; that is, 

S = G - B (A2) 

thus 

P = r <G - b) . 2 - g (as) 

If it is assumed that X = % and Y = §• 


then 


p = X - Y 


(A4) 


The standard deviation of p, cr(p), in terms of the standard deviations of X and Y is 

a(p) = [o2(X) + ct2(y3 1//2 (A 5) 


From the definitions of X and Y, their standard deviations can be given in terms of the 
standard deviations of K, G, and B; that is, 


similarly, 



(A6) 


(A7) 


Substitution of equations (A6) and (A7) into equation (A5) gives 


ct(p) = 


ct2(G) G 2 ct2(K) ot 2(B) B 2 cr2( K ) 

-f + JT + 


-.1/2 


K 


K 


K 




(A8) 




» 
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APPENDIX 


Division of equation (A8) by p yields the relative standard deviation or uncertainty in 
density 

g(p) _ 1 
P P 


a2(G) . G 2 a 2 (K) , ct 2 (B) , B 2 ct 2 (K) 

— 1 1 4* —— — -f- 


K* 


K 




K 


i/ * 


(A9) 


Substitution of equation (Al) into equation (A9) gives 


g(p) 

P 


K 

a 2 (G) 

s 

K 2 


| G 2 a 2 (K) t q2(B) f B 2 a 2 (K) 


-.1/2 




Hr 




or 


£(P) 

P 


q2(G) a2(B) q2(K)/G 2 + B 2 \ 

S 2 + S 2 + K 2 \ S 2 )\ 


(A10) 


Equation (A10) expresses the uncertainty in density in terms of the standard deviations of 
gross and background counting rates with respect to the signal and the standard deviation 
of the proportionality constant K modified by the term (G 2 + B 2 )/S 2 . 


The standard deviation of the gross or total counting rate is given by 



(All) 


where tg is the measurement time. This time is generally established by a trade-off 
between long time for minimum uncertainty in G (small cr(G)) and short time for mini- 
mum density change with altitude during the measurement period. 

The background counting rate can be assumed to consist of several factors: 

B = B d + B g + B n 


The standard deviation of B 


is given by 

q(B) = 



(A12) 


where tg is the time for background measurement. It is expected that ample measure- 
ment time will be available prior to atmospheric entry to determine the background 
counting rate with sufficient precision so that the uncertainty cr(B) will be negligible. 

The proportionality constant K can be stated in terms of a number of factors as 

K = f 1 f 2 f 3 t 4 R 


The calibration factor R is the counting rate per unit density of the sensor which is 
measured during sphere test calibration. The factors fj, f 2 , f3> and f 4 express 
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differences in the various parameters affecting sensor response between the time of 
measurement during atmospheric entry and the time of sensor calibration. For example, 

fl = the source emission rate at the measurement time compared with the emission 
rate during calibration 

f 2 = the amount of attenuation in the gamma- ray path at the measurement time com- 
pared with the attenuation during calibration 


f 3 = detection- system response at measurement time compared with the response 
during calibration 

f 4 = the atmospheric scattering efficiency, or cross section per unit density at 
measurement time, compared with the scattering efficiency during calibration. 

The relative standard deviation of K is then given by 


g(K) 

K 


a 2 (fl) | q2(f 2 ) [ pg^) | p2 ( f 4 ) 
f l 2 + f 2 2 + f 3 2 + f 4 2 



(A13) 


The factor fj is a function of the initial source strength, the half-lives of the pri- 
mary and any contaminant gamma-ray energies, and the source self-absorption. If the 
source to be used for flight can be used for calibration, then only the half-lives of the 
gamma-ray energies would need to be known to minimize o(f-jJ. The half-lives can be 
determined by taking a small sample of the main source and monitoring its output with a 
stable detection system during the flight mission. However, if different sources must be 
used for calibration and measurement, then the ratio of the gamma- ray outputs of the two 
sources must be precisely known. This ratio is determined from mass analysis of the 
flight and calibration sources which have been fabricated to have the same self- absorption 
characteristics and to contain radioactive material from the same production batch. In 
addition, direct-transmission counting- rate measurements of the gamma-ray output of 
each source must be made. 

The uncertainty in the amount of attenuation in the gamma- ray path is primarily 
a function of the shock layer formed outside of the spacecraft during entry. Even under 
extreme conditions of a thick, high-density shock layer, the attenuation in the shock layer 
is estimated to be less than 0.1 percent for a 100-keV gamma ray. Therefore, the atten- 
uation in the shock layer can be neglected. 

The factor fg is a function of the system parameters, gain, detector efficiency, and 
so on. Changes in these parameters such as gain shifts due to temperature or aging of the 
scintillation crystal will increase the uncertainty. However, as discussed in the text, the 
signal conditioning system will incorporate automatic gain control (AGC) circuitry which 
will compensate for changes of this type. Therefore, the uncertainty cr^fg) should be 
very small. 
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The factor is a function of the atmospheric scattering efficiency or cross sec- 
tion. As discussed in the text, this parameter varies only slightly with atmospheric com- 
position; therefore, the uncertainty 0^4) should be very small even for atmospheres 
containing appreciable argon. 

The parameter R, the counting rate per unit density, is determined by calibrating 
the sensor in a large vacuum sphere. The calibration experiment must be carefully 
designed to minimize extraneous effects such as scattering from the walls of the sphere 
so that R can be precisely determined. 
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TABLE 1.- MARS MODEL ATMOSPHERES 


[From reference 6) 


Property 

Symbol 

Dimension 

VM-1 

VM-2 

VM-3 

VM-4 

VM-5 

VM-6 

VM-7 

VM-8 

VM-9 

VM-10 

Surface Pressure 

p 0 

mb 

7.0 

7.0 

10.0 

10.0 

14.0 

14.0 

5.0 

5.0 

20.0 

20.0 



lb/ft 2 

14.6 

14.6 

20.9 

20.0 

29.2 

29.2 

10.4 

10.4 

41.7 

41.7 

Surface Density 

p o 

(gm/cm3)105 

0.955 

1.85 

1.365 

2.57 

1.91 

3.08 

0.68 

1.32 

2.73 

3.83 



(slugs/ft3)10^ 

1.85 

3.59 

2.65 

4.98 

3.7 

5.97 

1.32 

2.56 

5.30 

7.44 

1 Surface Temperature 

T o 

°K 

275 

200 

275 

200 

275 

200 

275 

200 

275 

200 



OR 

495 

360 

495 

360 

495 

360 

495 

360 

495 

360 

Stratospheric Temperature 

T s 

°K 

200 

100 

200 

100 

200 

100 

200 

100 

200 

100 



°R 

360 

180 

360 

180 

360 

180 

360 

180 

360 

180 

Acceleration of Gravity at Surface 

g 

cm/sec 2 

375 

375 

375 

375 

375 

375 

375 

375 

375 

375 



ft/sec 2 

12.3 

12.3 

12.3 

12.3 

12.3 

12.3 

12.3 

12.3 

12.3 

12.3 

Composition (percent) 













CO 2 (by mass) 



28.2 

100.0 

28.2 

70.0 

28.2 

35.7 

28.2 

100.0 

28.2 

13.0 

CO 2 by volume) 



20.0 

100.0 

20.0 

68.0 

20.0 

29.4 

20.0 

100.0 

20.0 

9.5 

N 2 (by mass) 



71.8 

0.0 

71.8 

0.0 

71.8 

28.6 

71.8 

0.0 

71.8 

62.0 

n 2 (by volume) 



80.0 

0.0 

80.0 

0.0 

80.0 

32.2 

80.0 

0.0 

80.0 

70.5 

Ar (by mass) 



0.0 

0.0 

0.0 

30.0 

0.0 

35.7 

0.0 

0.0 

0.0 

25.0 

Ar (by volume) 



0.0 

0.0 

0.0 

32.0 

0.0 

38.4 

0.0 

0.0 

0.0 

20.0 

Molecular Weight 

M 

mol" 1 

31.2 

44.0 

31.2 

42.7 

31.2 

36.6 

31.2 

44.0 

31.2 

31.9 

Specific Heat of Mixture 

Cp 

cal/gm°C 

0.230 

0.166 

0.230 

0.1530 

0.23 

0.174 

0.230 

0.166 

0.230 

0.207 

Specific Heat Ratio 

Of 


1.38 

1.37 

1.38 

1.43 

1.38 

1.45 

1.38 

1.37 

1.38 

1.41 

Adiabatic Lapse Rate 

r 

°K/km 

-3.88 

-5.39 

-3.88 

-5.85 

-3.88 

-5.14 

-3.88 

-5.39 

-3.88 

-4.33 



°R/1000 ft 

-2.13 

-2.96 

-2.13 

-3.21 

-2.13 

-2.82 

-2.13 

-2.96 

-2.13 

-2.38 

Tropopause Altitude 

hrp 

km 

19.3 

18.6 

19.3 

17.1 

19.3 

19.4 

19.3 

18.6 

19.3 

23.1 



kilo ft 

63.3 

61.0 

63.3 

56.1 

63.3 

63.6 

63.3 

61.0 

63.3 

75.8 

Inverse Scale Height (stratosphere) 

0 

km" 1 

0.0705 

0.199 

0.070 

0.193 

0.0705 

0.1655 

0.0705 

0.199 

0.0705 

0.145 



ft" 1 x 105 

2.15 

6.07 

2.15 

5.89 

2.15 

5.05 

2.15 

6.07 

2.15 

4.41 

Free Stream Continuous Surface 

V 

ft/sec 

186.0 

186.0 

156.0 

156.0 

132.0 

132.0 

220.0 

220.0 

110.0 

110.0 

Wind Speed 













Maximum Surface Wind Speed 

v 

max 

ft/sec 

380.0 

380.0 

310.0 

310.0 

270.0 

270.0 

450.0 

450.0 

225.0 

225.0 

Design Gust Speed 

v g 

ft/sec 

200 

200 

150.0 

150.0 

150.0 

150.0 

200.0 

200.0 

100.0 

100.0 



TABLE 2.- CALCULATED 1-SIGMA UNCERTAINTY IN 
GAMMA-RAY-BACKSCATTER MEASUREMENT OF 
ATMOSPHERIC DENSITY FOR MARS MODEL 
ATMOSPHERES VM-8 AND VM-9 


Altitude 

Uncertainty in density, percent 

ft 

km 

VM-8 

VM-9 

0 

0 

3.2 

3.1 

5 000 

1.5 

3.2 

3.1 

10 000 

3.1 

3.2 

3.1 

15 000 

4.6 

3.3 

3.1 

20 000 

6.1 

3.3 

3.1 

25 000 

7.6 

3.4 

3.1 

30 000 

9.2 

3.4 

3.1 

40 000 

12.2 

3.6 

3.2 

60 000 

18.3 

4.3 

3.2 

80 000 

24.4 

7.0 

3.3 

100 000 

30.5 

15.3 

3.5 

120 000 

36.6 

41.5 

3.7 

140 000 

42.7 

133.0 

4.1 

160 000 

48.8 

435.0 

4.6 

180 000 

54.9 


5.3 

200 000 

61.0 


6.5 

220 000 

67.1 


8.0 

240 000 

73.2 


10.3 

260 000 

79.3 


14.0 

280 000 

85.4 


19.1 

300 000 

91.5 


26.9 

320 000 

97.6 


38.3 

340 000 

103.7 


55.7 

360 000 

109.8 


81.9 

380 000 

115.9 


121.0 

400 000 

122.0 


181.0 
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Figure 2.- Relative probability of occurrence of 
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Figure 3.- Payload of Nike-Apache backscatter density sensor. 
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Detected counting rate (background subtracted) and density as a function of altitude for an earth-atmosphere flight test 

(data from ref. 5). 
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Figure 5.- Density as a function of altitude for extremes of Mars model ati 
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Figure 7.- Calculated 1-sigma uncertainty in density measurement as a function of altitude for extremes of model atmospheres, VM-8 and VM-9. 
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